Coherent versus Incoherent Light Scattering from a Quantum Dot 
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We analyze the light scattered by a single InAs quantum dot interacting with a resonant 
continuous- wave laser. High resolution spectra reveal clear distinctions between coherent and in- 
coherent scattering, with the laser intensity spanning over four orders of magnitude. We find that 
the fraction of coherently scattered photons can approach unity under sufficiently weak or detuned 
excitation, ruling out pure dephasing as a relevant decoherence mechanism. We show how spectral 
diffusion shapes spectra, correlation functions, and phase-coherence, concealing the ideal radiative ly- 
broadened two-level system described by Mollow. 

PACS numbers: 78.47.-p, 78.67.Hc 



Like an isolated atom or ion, a semiconductor quan- 
tum dot (QD) "artificial atom" scatters monochromatic 
laser light incident upon it. For an ideal two- level system 
broadened by radiative decay at a rate ft, the spectral and 
temporal properties of the scattered photons are deter- 
mined solely by the laser detuning from resonance, Acj, 
and by the Rabi frequency, Q. In seminal work, Mollow 
showed that when Q <C ft, most of the light is scattered 
coherently, i.e. elastically, whereas it is otherwise dom- 
inated by resonance fluorescence^. Moreover, when the 
scattered light originates from a single two-level system 
it exhibits photon anti-bunching, for any value of Q and 
AuP. 

Probing and controlling quantum phase coherence is 
at the core of quantum information science, and semi- 
conductor QDs are well suited for investigating quan- 
tum coherence in solids at optical frequencies^. In this 
context, resonant light scattering in QDs has been of 
special interest, and has emerged as a promising re- 
source for the generation of highly ideal single pho- 
ton states^. Milestone demonstrations include oscilla- 
tory field correlations^, Mollow triplets 6 8 , photon anti- 
bunching and cascaded photon emission 9 . Nevertheless, 
the theoretical modelling of these experimental obser- 
vations has necessitated the inclusion of a phenomeno- 
logical dephasing time, T 2 , in general smaller than 2Ti, 
where T\ = 1/ft is the radiative lifetime^ 9 -] Moreover, 
the coherent scattering — as opposed to incoherent scat- 
tering that causes resonance fluorescence — has remained 
largely unexplored. A more detailed experimental inves- 
tigation clarifying the role of decoherence in the context 
of resonant light scattering is all the more timely given 
recent proposals of using coherent scattering for the gen- 
eration of indistinguishable photons^SEil theory predicts 
that dephasing reduces the fraction of coherently scat- 
tered photons from unity, regardless of Q and Auj. Here 
we combine high resolution spectroscopy (35 MHz), pho- 
ton correlations and phase-coherence measurements, to 
obtain clear distinctions between coherent and incoherent 
scattering. We show that the previously observed non- 



ideal linewidth and correlation functions are due to the 
fluctuation of the QD resonance frequency on a time scale 
longer than Ti , resulting in an apparent broadening of the 
lineshape. One the other hand, at the relevant time scale 
of quantum evolution, the system is close to an ideal two- 
level system with negligible pure dephasing, as described 
by Mollow^, and as expected fr om four- wave mixing mea- 
surements on QD ensemble d 12 * 13 *. We extend Mollow's 
theory to include spectral diffusion in the form of inho- 
mogeneous broadening and provide a complete picture of 
resonant light scattering in semiconductor QDs. 

We probe QDs grown by molecular beam epitaxy at the 
center of a planar optical microcavity. QDs of this type 
have been investigated extensively due to their atom-like 
spectrsP and long dephasing times^. The QD sample 
was grown using a solid source VEECO Gen-II molec- 
ular beam epitaxy (MBE) system on a semi-insulating 
GaAs (100) substrate with a 12 pairs top and 20 pairs 
bottom Alo.9Gao.1As/GaAs distributed Bragg reflector. 
The substrate rotation was stopped during MBE growth 
of InAs layers to obtain a QD density varying uniformly 
from 10 9 /cm 2 to 10 8 /cm 2 ^'. The dominant vertical cav- 
ity mode is centered around A^925 nm. The sample was 
maintained at a temperature of 3.8 K in a closed-cycle 
cryostat and the QD emission was collected by an in situ 
high numerical aperture aspheric lens [Fig. 1(a)]. To ex- 
cite a QD resonantly and discriminate the resonant scat- 
tering signal from stray laser light .an orthogonal excita- 
tion/detection geometry was usecP. The detected light 
was coupled into a single mode fiber without prior filter- 
ing or cross-polarized excitation/detection. For coarse 
spectral analysis of the emitted light a grating spectrom- 
eter with a cooled charge-coupled device camera was em- 
ployed. For high resolution spectral measurements we 
used a scanning Fabry-Perot interferometer with a free 
spectral range of 4.2 GHz and finesse of 120 in conjunc- 
tion with a single photon counting detector. The excita- 
tion source was a diode laser continuously tunable over a 
range of several GHz. We have probed a number of single 
QDs and found very similar behavior, although the data 
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FIG. 1: (Color online) (a) Schematic of experimental setup, 
(b) QD emission spectrum, recorded with a grating spectrom- 
eter under resonant laser excitation as is (black solid trace), 
20 x magnified (dashed gray trace) and on a logarithmic scale 
(solid red trace) to visualize the phonon broadband, (c) QD 
excitation spectrum, (d) Spectrally integrated intensity of 
scattered light as a function of time showing flickering (dark 
blue trace) that is inhibited when an additional weak auxil- 
iary laser is added (light green trace). The black trace was 
recorded with the auxiliary laser only. 



presented here is from one specific QD. 

Figure 1(b) shows the power spectrum of this QD 
under resonant laser excitation, recorded with a spec- 
trometer with ^10 GHz resolution. Although the main 
emission line is not resolved in this measurement, a 
broadband emission around it can be clearly identified. 
As is well-known from experimental photoluminescence 
(PL) measurements^, four-wave mixing studies^^, and 
from theoretical calculations of resonance fluorescence 
spectra^, this broad emission originates from fast scat- 
tering processes with acoustic phonons. Although at 3.8 
K [Fig. 1(b)] this phonon scattering "pedestal" is highly 
asymmetric it becomes mor e sym metric and prominent 
with increasing temperatur e 1 12 * 15 H£J In Fig. 1(b) about 
5% of light is emitted into this broad band, thus at liq- 
uid He temperature as much as ^95% of light may be 
scattered coherently, irrespective of Q and Auj. 

Figure 1(c) shows the QD excitation spectrum at an 
excitation intensity below the saturation intensity, ob- 
tained by scanning the laser across the QD resonance 
frequency while collecting the scattered light. The full 
width at half maximum (FWHM) of the resonance is 
about 0.7 GHz. Although in four- wave mixing stud- 
ies the zero-phonon linewidth has been shown to be 
limited by radiative decay to about 170 MH^PHH, PL 
linewidths and/or resonantly measured single QD exci- 
tation linewidths are usually found to lie in the range 
of 500 MHz to several GHz ^ 11 ' 17 ' 19 ^ . The additional 
"apparent" broadening is due to spectral diffusion, a pro- 



cess by which the QD transition frequency is randomly 
shifted during the measurement 17 . This shift is thought 
to originate from a fluctuating charge environment of the 
QD that occurs on a time scale that is long compared to 
the radiative decay process. In PL measurements spec- 
tral diffusion is manifested in the form of long-time pho- 
ton correlations and blinking 22 . In resonant scattering 
experiments the existence of charge fluctuations is evi- 
denced by nickerin^pU and by the fact that many res- 
onances identified in PL as neutral exciton transitions 
generate little or no resonance fluorescence unless a weak 
background lighting is presenP°4 It has been proposed 
that a Coulomb blockade effect is at the origin of this res- 
onant emission quenching 20 . Flickering of the resonant 
scattering signal is directly observed when no external 
non-resonant lighting is present, as seen in Fig. 1(d), 
consistent with spectral diffusion 24 . The flickering can 
be inhibited on the measurement timescale (~10 ms) by 
an auxiliary non-resonant light source (here a laser at 
A=660 nm) that is so weak that by itself generates neg- 
ligible PL [Fig. 1(d)]. We speculate that the effect of 
this auxiliary laser is to change the timescale of the flick- 
ering by neutralization of surrounding charges. Spectral 
diffusion is also commonly observed for QDs in nanos- 
tructures, where it is associated with the proximity of 
etched surfaces 2 ^. 

In order to distinguish experimentally between coher- 
ent and incoherent scattering it is necessary to analyze 
the scattered light spectrally with a resolution better 
than n. For a two- level system with natural resonance 
frequency Co>o, exposed to a monochromatic field of fre- 
quency lj, the power spectrum of the scattered light takes 
the analytic form, given by MolWD 



g(v,u),LJo) = 27r|a 00 | 2 <5(i;-w)+n 00 Kf2' 



\Mv-"W 
(1) 

where v is the emission frequency, and the only source of 
broadening is a decay of the upper state to the lower state 
at a rate k due to spontaneous emission. The first term 
corresponds to photons scattered coherently while the 
second term describes the resonance fluorescence. The 
laser detuning Auj = uj — ujo enters through the poly- 
nomial function /, the steady- state population inversion 
noo, and the steady-state quantum mechanical expecta- 
tion value of the two- level coherence, a^, given inland in 
Appendix A. Note that the function g(v,uj,uJo) is always 
symmetric around v = uj (not v = cjq), unlik e spectra of 
dressed states populated non-resonant ljffHUD 

Figure 2(a) shows maps of the scattered light as a func- 
tion of emission and excitation frequency, for various val- 
ues of the Rabi frequency. The latter was varied by vary- 
ing the laser intensity. As expected from Eq. (1), when 
Q <C n [leftmost map in Fig. 2(a)], the emission reso- 
nance linewidth is much less than n. It is limited here by 
the resolution of our scanning Fabry-Perot interferome- 
ter (35 MHz) but is expected to be as narrow as the laser 
linewidth (~1 MHz). With increasing Q the resonance 
fluorescence eventually dominates in the form of the Mol- 
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FIG. 2: (Color online) (a) Maps of scattered light intensity as a function of detection frequency (abscissas) and excitation 
frequency (ordinates) relative to the QD transition frequency, for three different values of the Rabi frequency. Each panel was 
recorded in 200 s. At Q/27r=0.13 GHz the detector count rate was 3xl0 5 s _1 at the input of the Fabry-Perot interferometer. 
The dashed lines indicate the location of the eigenfrequencies of the coupled laser/QD system. The inset shows the line section 
v — uj together with the corresponding theoretical curve, obtained from Eq. (2). The faint diagonal lines parallel to the v — uj 
section are due to the residual transmission at high-order modes of the Fabry Perot interferometer and satellite modes of the 
laser, (b) Theoretical maps [Eq. (2)] corresponding to a radiatively-broadened two-level system subject to spectral diffusion, 
(c) Theoretical maps [Eq. (1)] representing an ideal, radiatively-broadened two-level system, (d) Theoretical maps [Eq. (All)] 
representing a two- level system subject to radiative decay, and to pure depahsing at a rate 7. 7 is chosen so as to obtain the 
excitation linewidth of Fig. 1(c). 



low triplet, composed of three peaks with a FWHM that is of the order of k [middle and rightmost map in Fig. 
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2(a)]. 

In order to identify the dominant effects on the QD 
two-level system due to its solid-state environment, we 
make a side-by-side comparison of our experimental data 
of Fig. 2(a) with spectral maps obtained from three 
different theoretical models, displayed in Fig. 2(b-d). 
In all cases, we assume a radiative decay, ft/27r=180 
MHz, which is consistent with typical radiative lifetimes 
[t=1/(2ttx180 MHz)=0.9 ns] for InAs QDs. We include 
spectral diffusion by integrating Eq. (1) over all possible 
(random) detunings to obtain 



where we assume a Gaussian distribution of QD reso- 
nance frequencies due to spectral diffusion with FWHM 
s « 2.355<j. While spectral diffusion may manifest at dif- 
ferent timescales^, we assume here that this timescale is 
long compared to the radiative decay time. We also ne- 
glect any influence of dark excitonic stateiPH To account 
for finite apparatus resolution we replace the delta func- 
tion in Eq. (1) by a normalized Lorentzian with FWHM 
equal to 35 MHz. Fig. 3(b) shows spectral maps ob- 
tained using Eq. (2) with ,s/27r=0.7 GHz, chosen to co- 
incide with the measured excitation linewidth of Fig. 1(c) 
which is also typical for InAs QDiP^^ ^ 21 1 23 1 . In com- 
parison, Fig. 2(c) shows spectral maps obtained using 
Eq. (1) clearly not agreeing as well with the experimen- 
tal data of Fig. 2(a) and highlighting the subtle effects 
of the spectral diffusion process. In Fig. 2(d) we further 
show the case when pure dephasing is present instead of 
spectral diffusion, as was assumed in previous studied^. 
The most general theoretical expressions for the power 
spectrum for the latter case are derived in Appendix A. 
The pure dephasing rate, 7, was chosen to yield the ob- 
served excitation linewidth [Fig. 1(c)]. As is evident from 
this comparison, pure dephasing does not play a major 
role at the measurement temperature since it would give 
rise to rather different observations. In particular, the 
absence of broad emission superimposed to the leftmost 
spectral map in Fig. 2(a) is an unequivocal indication 
that dephasing is not a relevant broadening process here. 

For a quantitative analysis we recorded a series of spec- 
tra at exact resonance (Acj = 0) with sufficiently high 
signal to noise ratio for precise comparison with theory 
(Fig. 3). The same information is displayed on both 
linear (left) and logarithmic (right) ordinate scales. The 
theoretical traces were obtained using Eq. (2) with the 
same parameters as in Fig. 2(b), and only a common 
scale factor was permitted for all traces. Short-dashed 
(long-dashed) red lines correspond to the total (incoher- 
ent) scattered light intensity and the light-red shaded 
area indicates the coherently scattered light. As is seen 
in the figure, there is excellent agreement between the 
experimental data and the theoretical expression of Eq. 
(2) evaluated numerically, considering that the excitation 
laser intensity spans more than four orders of magnitude. 
Since at the largest Q in Fig. 3 the central peak continues 



to grow, some amount of light is at that point originat- 
ing from other scatterers, perhaps neighboring detuned 
QDs. Nonetheless, when Q/2tt=0A9 GHz, we can esti- 
mate that more than 90 % of the detected light originates 
from scattering off the QD probed. 
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FIG. 3: (Color online) Power spectrum of light scattered by 
the QD at exact resonance (Acj = 0) represented on a linear 
(left) and logarithmic (right) ordinate scale, for a range of 
Rabi frequencies. Each spectrum was recorded in 60 s. The 
theoretical curve (red dashed line) was obtained by numerical 
evaluation of Eq. (2). The light red shaded area corresponds 
to coherently scattered light which dominates whenever Q < 
k. If pure dephasing had a significant role in the scattering 
process, a broader feature would always be visible, even when 
Q < k. 
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Photon statistics of the scattered light, shown in Fig. 
4(a) below saturation (left panel), near saturation (mid- 
dle panel) and above saturation (right panel), corrob- 
orate this via the observation of photon antibunching. 
Photon correlation measurements were performed with a 
Hanbury-Brown and Twiss setupP^that uses a 2x2 fiber 
coupler and two fiber coupled avalanche photon detec- 
tors. The detectors are nominally identical to those used 
iiP and thus an identical instrument response function 
was assumed. We show the raw unprocessed normalized 
data as well as the theoretical normalized second order 
correlation function in the presence of spectral diffusion, 

(2) 

9sr>(t)> convolved with the instrument response function 
(IRF) 6 . It is the IRF, not stray background light, that 

(2) 

causes the dip to rise above zero near t=0. 9sr>(f) is com- 
puted identically to Eq. (2), replacing g(y, uj' ) under 
the integral with g( 2 \t,u,u , Q)n 00 /n 00 {uj = ujq), where 
g( 2 \t,oj,oj' ) is the normalized the correlation function 
for the ideal two-level system. Above saturation Rabi 
oscillations, the time-domain analogue of the side bands 
of Fig. 2, are seen 6 . 

Spectral diffusion significantly alters the spectral and 
temporal characteristics of the scattering process com- 
pared to the ideal two- level system. Figure 4(b) shows 
the total normalized scattered light intensity as a func- 
tion of Rabi frequency for the parameters used in Fig. 3. 
using Eq. (2) (solid red trace) and Eq. (1) (solid blue 
trace). Correspondingly, the dashed lines of the same 
color represent the fraction of the intensity of coherently 
scattered light alone. The same are plotted in Fig. 4(c) 
but under significant laser detuning. A signature effect of 
the spectral diffusion is to cause an increase in the frac- 
tion of photons scattered coherently, as would any source 
of inhomogeneous broadening. At large values of f], more 
light is actually scattered coherently with the laser off- 
resonance than with the laser at exact resonance. This 
is seen in Fig. 4(c) but is also directly visible in the bot- 
tom map of Fig. 2 where a minimum occurs along the 
line v — uj at resonance (inset). 

Lastly we examine the phase coherence between the 
laser and the scattered light (below saturation) by com- 
bining the scattered light with a local oscillator (LOSC) 
signal at a beam splitter. We measure the fringe contrast 
obtained when varying the LOSC phase with a piezoelec- 
tric actuator (PZT) [Fig. 4(d)]. The fringe contrast we 
obtain here is ~40% as seen in Fig. 4(e) which can be 
understood by spectral diffusion that causes large fluc- 
tuations in photon flux at the beam splitter. In theory 
we expect spectral diffusion to reduce the fringe contrast 
(visibility) by a factor of order n/ s ^0.25. 

In summary, we have carried out high-resolution mea- 
surements of the light scattered coherently and incoher- 
ently by a single In As QD, revealing in great detail how 
the scattering process evolves over more than four orders 
of magnitude of excitation laser intensity. The simple 
inclusion into Mollow's theory of spectral diffusion as a 
source of inhomogeneous broadening does faithfully re- 
produce both our observed spectra but also photon cor- 
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FIG. 4: (Color online) (a) Photon statistics of scattered light 
for three different Rabi frequencies. When Q > k, Rabi oscil- 
lations appear, (b) Plot of the total scattered light intensity 
(solid red trace) and the fraction of coherently scattered light 
intensity (dashed red trace), using the same parameters as 
those used in Fig. 3. For comparison the same are plotted 
in blue when not including spectral diffusion, (c) Same as in 
(b) but with a laser detuned by Aoj/2ty = 1 GHz. (d) Mea- 
surement of mutual phase coherence between the coherently 
scattered light and a local oscillator (LOSC) by interferome- 
try. (e) Intensity of light at the output of the beam splitter 
in (d) as a function of LOSC phase. 



relation and phase-coherence measurements. The insight 
that an apparent broadening rather than pure dephasing 
predominantly affects the scattering process has impor- 
tant implications for future use of QDs, for example in 
quantum repeaters^. There we expect spectral diffusion 
to reduce the photon flux without however limiting two- 
photon indistinguishability^. The temporal flickering of 
the resonant scattering signal associated with the spec- 
tral diffusion may be reduced in future improved struc- 
tures. For example it may be possible to control the 
timescale for this flickering by placing QDs in charge- 
tunable devices. Overall our work offers a complete pic- 
ture of resonant light scattering relevant to a wide variety 
of solid-state nanostructures. 

The authors acknowledge financial support from the 
National Science Foundation (NSF DMR-0906025 and 
CMMI-0928664) the National Natural Science Founda- 
tion of China (Grant No. 90921015). 
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Appendix A: Power spectrum in the presence of 
pure dephasing and detuning 

Using the same notation as in RefpJ we provide an 
expression for the power spectrum of the light scattered 
by a two- level system subject to an incident monochro- 
matic laser, when including an additional off-diagonal de- 
cay (pure dephasing) rate 7. Unlike in Ref.^we further 
assume that the detuning is not zero, i.e. Auo ^ 0. 

The optical Bloch equations in the rotating-wave ap- 
proximation read 



d 
dt 



R(t) = M • R(t) (Al) 



in which 



M = 



and 



-k -in/2 in/2 

-in/2 + in/2 

in/2 ' -tf/2-7 -iAuj -in/2 



in/2 



-in/2 



R(t) 



/ n(t) 
a(t) 
a*(i) 

\ m(t) 





(A2) 



(A3) 



Here n(t)=Tv{p(i)a) a}, a(t)=Tr{p(t)a}, 

a*(t)=Tr{p(t)aJ}, and m(t)=Tr{p(t)aa^}, where a, 
a\ and p are the lowering, raising, and density op- 
erators, respectively. The diagonal and off-diagonal 
damping terms ^- = k and = 7 + ^/2, respectively, 
have been included using the usual master equation. 
Equation I Al I has the steady state solution 



in 



k + 27 + 2iAcu 



4 Auj 2 + (k + 2 7 ) (k + 2 7 + 2ft 2 /«) /4 



and 



o 2 



K + 27 



4k Aw 2 + (k + 27) (k + 27 + 2fi 2 //t)/4 



(A4) 



(A5) 



Following Mollow^, the power spectrum is obtained as 
the Fourier transform of the first-order correlation func- 
tion, g(r,t) = {a)(t)a(t + r)), in steady-state (t — » 00), 
as 

/oo 
ff (r)e il ' T dr = 2Re(5(-i^)) (A6) 
-CO 



where is the angular frequency of the scattered light, 
and g(s) is the Laplace transform of g(r). In order to 
calculate g(r) we make use of the quantum regression 
theorem^ which states that for an operator O whose 
expectation value is known to evolve from time t to time 
t + r as 



<0(t + T))=$>ito<0i(t)> 



(A7) 



the two-time correlation function (Oi(t)0(t + r)Ok{t)) 
can be calculated as a function of single time expectation 
values as follows: 

<Oi(t)0(t + T)O fc (t)> = Y.^iOmO^Okit)) 

(A8) 

^Given that the solution of equation Al can be written as 
K(t + r) = e Mr • R(t) (A9) 
we obtain 



5(T) = (e MT |2,2noo+e MT | 2 ,4aye 



Mr I * \ —iujT—S F pr/2 



(A10) 



To account for the limited resolution, 5fp, of our scan- 
ning Fabry- Perot interferometer, we have multiplied the 



A6 



correlation function by e~ 6FpT / 2 (r > 0). Computing 
requires finding the Laplace transform of |A1Q[ which in 
turn involves computing the matrix (Is — M) _1 , where I 
is the 4x4 identity matrix. We finally obtain 

g{y) = 2Re((I(-w/ + iuo + 5 F p/2) - M) -1 ^ 



(!(-«/ + zcj + Sfp/2) - M)- 1 ^,^; 



(All) 



which for the special case 7 = and Spp = reduces to 
the expression 



)(v) = 2Tr\a O0 \ 2 S(iy — u) + UootiVf 



{v - uj) 2 + Sl 2 /2 + k 



\f(i(v-co)W 

(A12) 

given by Mollow [Eq. (1)], where f(s) = s 3 + 2ks 2 + 
(n 2 + (Aw) 2 + (5/4)k 2 )s + k(\<D? + (Aw) 2 + \k 2 ). 
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